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Abstract. We describe new implementations of quantum error correction that are continuous in time, and thus described by 
continuous dynamical maps. We evaluate the performance of such schemes using numerical simulations, and comment on the 
effectiveness and applicability of continuous error correction for quantum computing. 



INTRODUCTION 

Error correction is an essential ingredient for quantum computing and communication. The combined use of error 
correction and fault tolerance techniques allows a quantum computation to succeed even in the presence of errors 
during the computation. Quantum error correction is possible because of the ability of quantum error correcting codes 
(QECC) to redundantly encode quantum information and permit its retrieval. The first of these codes were discovered 
by Shor [1] and Steane [2], and these discoveries sparked a slew of work in quantum error correcting codes, part of 
which has resulted in the elegant theory of stabilizer codes [3]. Stabilizer codes are the best known (class of) QECC 
to date and most error correction protocols involve their use ' . Error correction schemes that involve the use of QECC 
and the detection and correction operations that go along with them are refeiTed to as active error correction schemes. 

Figure 1 illustrates the standard blueprint for implementing active eiTor correction schemes. The procedure begins 
by encoding of the (unknown) qubit state to be protected using an error connecting code. Then after a period of 
time during which errors could occur, an eiTor detection operation is performed which, assuming the code used 
can handle the number and type of errors that has occuiTed, yields information about eiTor type and location. This 
information is used to quantum mechanically or classically condition the next block - the correction - which performs 
operations on the encoded block of qubits to reverse the effect of the eiTor. If the conditioning between the detect and 
correct blocks is quantum mechanical, then this is referred to as eiTor connection without measurement. If on the other 
hand, the conditioning is classical then it is referred to as error correction with measurement, and requires classical 
communication between the two blocks. 

The detection and correction operations are often repeated many times, and the whole sequence concludes with 
a decoding operation that retrieves the protected one qubit state from the encoded block. A feature of this blueprint 
we would like to emphasize is its discrete nature. The detection and correction operations occur at discrete points in 
time and they are usually assumed to be instantaneous and arbitrarily accurate. For this reason we will refer to such an 
implementation model as discrete error correction. In this paper we will present an alternative to this model where both 
the detection and correction blocks operate continuously and occur at a certain rate as opposed to instantaneously. This 
will lead to a description of error correction in terms of continuous dynamical maps, which we will solve numerically 
to evaluate error correcting capabilities. 

There have been a number of authors in the past who have considered continuous forms of eiTor correction. Paz and 
Zurek [5] examined the limit where the detection and correction operations are performed directly after each other 
continually. However, this is not quite continuous error coiTection in the sense of this paper because they still assumed 
both operations occur instantaneously. Mabuchi and Zoller [6], and Ahn et. al. [7] have examined continuous error 
correction for a specific subset of error models (detected errors). Mabuchi and Zoller also assume that the correction 



' For completeness, we mention that there are also other techniques that protect quantum information against errors (e.g. decoherence free 
subsystems [4]). 




FIGURE 1. Standard active error correction implementation blueprint. The communication between the detection and correction 
blocks will require classical signals for error correction with measurement. 

operations are instantaneous, and therefore treat a slightly different case from us. Finally, Ahn et. al in [8] examine the 
exact case we will cover in the section entitled continuous error correction by indirect feedback. However, the error 
correction scheme they describe involves intensive real-time computation, and is therefore considerably less practical 
than the scheme we detail. 



REVffiW OF CONCEPTS 
Quantum feedback control 

Making the detection and correction operations continuous places the process of error correction within the frame- 
work of quantum feedback control [9, 10, 11, 12, 13] - a field which examines the control of quantum systems using 
closed loop control methods. We will use the tools of this field to describe the error correction procedures below. 

We will not cover the field of quantum feedback control in any detail here, for that the reader is referred to [13, 12] 
and references therein. However, it is useful for what follows to make a distinction between the two types of quantum 
feedback control. This distinction arises essentially because the controller steering the quantum system of interest 
can itself be either quantum or classical. The two alternatives lead to different control loops and hence different 
mathematical descriptions. This distinction is also presented in [12] and we shall use the terminology used there. 

A direct feedback control loop has a quantum controller that is directly coupled to variables of interest in the 
quantum system to be controlled. The interaction between the system and controller in this case is unitary. In an 
indirect feedback loop, measurements are performed on a quantum system and the results of these measurements are 
used by a classical controller to form classical conditioning signals to control the system. 



The bit-flip code 

As mentioned earlier, stabilizer codes are the best known QECCs. We will not attempt a survey of the extensive 
theory of stabilizer codes (the interested reader is referred to [3, 14]), but will rather describe one of the simplest 
stabihzer codes - the bit-flip code - which will be used in the following sections to illustrate the continuous error 
correction schemes. 

The bit-flip code is essentially a classical parity check code. It protects against a logical flip error which takes one 
computational basis state to the other. The bit-flip code protects against this error by using the following repetition 
encoding: |0)^ = |000)p, |1)^ = |111)/), where the subscripts L and P stand for logical and physical, respectively. 
Therefore a general encoded qubit will have the form | i/a) = a |0)^ + j3 1 1)^ = a |000) + j3 1 1 1 1). We will refer to the 
encoded qubit states as the codewords, and the subspace they span as the codespace. 

This code can detect and correct one bit-flip. The detection operation involves measuring the operators ZZI and IZZ 
^, which are referred to as the error syndromes or stabilizer generators. Two things to note, both of which are properties 



^ We denote the Pauli Ox,oy, and Oz operators by X, Y, andZ, respectively, and suppress the tensor product sign. Therefore ZZI = Z®Z®1. 



TABLE 1. The three qubit bit-flip code. Note that each 
error results in a different sequence of stabihzer generator 
measurement results. (■)p represents the expectation value 
of • under the encoded three qubit state p. 
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FIGURE 2. Circuits for implementing the three qubit bit-flip code. The circuit on the left shows an implementation that uses 
measurement, and the circuit on the right shows one that uses no measurement. In both figures, the top three qubits form the 
encoded logical qubit and the bottom two are ancilla. Note that for the circuit on the right, to repeat the error correction procedure, 
the ancilla qubits must be replaced or reset to the |0) state at the end of each run (at the far right of the circuit). 



of general stabilizer codes, are that the stabilizer generators commute, and that the codewords are both eigenvalue one 
eigenstates of the stabilizer generators (or in other words, the codespace is stabilized by the generators). 

The four possible outcomes from measuring the two stabilizer generators tell us about the four possible error events. 
This is illustrated by table 1. Correcting errors using this code then simply amounts to applying a unitary to restore the 
encoded state back to its unperturbed value. The value of this unitary depends on the measurement results as table 1 
shows. 



CONTINUOUS ERROR CORRECTION BY INDIRECT FEEDBACK 

The indirect feedback model of error correction is a continuous version of discrete error correction with measurement. 
Full details on this model are in the paper [15], we will simply outline the results. The subfigure on the left of figure 
2 shows a circuit which implements the bit-flip code in a discrete manner with measurement. The top three qubits 
in the circuit represent the encoded qubit and the bottom two represent ancilla qubits used in the stabilizer generator 
measurements. The first part of the circuit, the CNOT gates and the projective measurements, detect the error by 
placing the results of the stabihzer generator measurements in the ancilla qubits and measuring in the computational 
basis. Then the results of these measurements are used to condition the unitary correcting gate, ^. To make this process 
continuous we must replace the CNOT gates and projective measurements by a weak, continuous measurement of the 
stabilizer generators, and the unitary correcting gate by Hamiltonian evolution which is conditioned on the continuous 
stabihzer measurement record. 

Note that for the three qubit code we need to measure two stabilizer generators continuously and simultaneously. 
This poses no fundamental problems because by definition the stabilizer generators are commuting observables. Also, 
because the codespace is stabilized by the generators, when there is no error, these measurements do not affect the 
encoded qubits at all. However, one problem that does arise is that the continuous measurements, because they are 
weak, are very noisy. This is just a manifestation of the tradeoff in quantum mechanics between the strength of a 
measurement (and hence the amount of back action on the system being measured) and the the amount of information 
gained from the measurement. As in classical feedback control, feedback conditioned on noisy measurement records 
is in general ineffective, and thus we cannot perform error correction by using the raw measurement signals. Instead 
we must insert a signal processing step between the measurement and the classical controller which smoothes the 
stabilizer generator measurement records. 



Therefore the steps involved in the error correcting scheme we propose are: 

1. Encode information in a stabilizer code suited to the errors of concern. 

2. Continuously perform weak measurements of the stabilizer generators, and smooth the measurement records. 

3. From these smoothed measurement records, form conditioning signals for feedback operators on each physical 
qubit. The mapping between measurement results and conditioning signals is determined by the stabilizer code. 

4. Apply feedback (correction) Hamiltonians to each physical qubit, where the strength of the Hamiltonians is given 
by the conditioning signals formed in the previous step. 

For the bit-flip code, weak measurements are performed of the observables ZZ/ and /ZZ. The feedback Hamiltonians 
are XII, IX I, and IIX, each one implementing a correction on one of the physical qubits in the encoding. The mapping 
between the results of the (smoothed) measurements and the feedback conditioning signals is essentially a translation 
of table 1. 

We use the language of open quantum systems [16, 17, 18] to model this process. The following stochastic master 
equation (SME) describes the evolution of the encoded three qubit system while undergoing random bit-flip errors, 
weak measurement of the stabilizer generators, and continuous feedback conditioned on the measurement results [15]. 

dpc{t) = y{S>[XII] + S![IXI] + &[IIX])p,{t)dt + K-(^[ZZ/] + 9[IZZ\)p,{t)dt + ^/K{.^[ZZI]dW\ (f) 

+ ^[IZZ]dW2{t))pc{t) - iX{Gi{t)[XII,pc{t)]+G2{t)[IXI,pc{t)]+Gi{t)[IIX,pc{t)])dt (1) 

where y is the error rate for each qubit, K is the measurement strength, and A is the maximum feedback strength. 
Gi{t) are the feedback conditioning signals formed from the smoothed measurement records, the quantities dWi{t) are 
Weiner increments [16], and and ^ are the superoperators 

&[A]p =ApA^ -^{A^Ap+pA^A), ^[A]p = Ap + pA^ - p XT[Ap + pA^] (2) 

for any operator A. Also note that we set Tz = 1 throughout this paper. 

Eq. (1) is a non-Markovian, non-hnear SME, and therefore the only way to solve it is through numerical simulation. 
We did precisely this (for details see [15]), to evaluate the efficacy of this error correction scheme. The results of these 
simulations are summarized by the left subfigureof figure 3. We used fidelity with the initial state, F(f ) = (i//o| p(f) | V^o) 
as a figure of merit to evaluate the error correction scheme. The figure shows fidehty versus time curves for several 
values of error rate (y). Each plot also shows the fideUty curve for one qubit in the absence of error correction. 
A comparison of these two curves shows that the fidelity is preserved for a longer period of time by the error 
correction scheme for small enough error rates. Furthermore, for small error rates (7 < 0.3) the solid curve shows 
a vast improvement over the exponential decay in the absence of error correction. However, we see that as the error 
rate increases, the fidelity decay even in the presence of error correction behaves exponentially, and the two curves 
look very similar. When the error rate gets large enough, the error correcting scheme becomes unable to handle the 
errors and becomes ineffective. 



CONTINUOUS ERROR CORRECTION BY DIRECT FEEDBACK 

The direct feedback model of error correction is a continuous version of discrete error correction without measurement. 
To illustrate this model we will again describe the continuous analogue of a circuit that implements the bit-flip code: 
the subfigure on the right of figure 2. In this circuit the results of the stabihzer generators are again placed in the ancilla 
qubits, however, instead of measurements followed by classical conditioning of the correction operation, the correction 
is done by direct coupling between the ancilla and the encoded qubits (via Toffoli gates in this case). It is important to 
note that the ancilla qubits must be reset to the |0) state after each run of the circuit. This is essentially a consequence 
of the no-deleting theorem [19] for quantum information and because the error is moved from the encoded qubits into 
the ancilla array. 

The process of making this circuit continuous is much simpler than in the indirect feedback case. We simply have 
to model the detection and correction unitary gates as one continuous Hamiltonian evolution. However in addition to 
this, we must also replace the ancilla reset procedure with a continuous map. The natural candidate is a cooling process 
that continuously damps each ancilla independently. Both these processes, along the with random bit-flip errors, can 




FIGURE 3. Fidelity curves for the two error correcting sctiemes. The figure on the left shows curves for several error rates 
for the indirect feedback model. The solid curves are the fidelities of the encoded qubit (three qubit code) with parameters: 
K = 150 Hz,?i = 150 Hz, initial state Iyo) = |000)). The figure on the right shows curves for several Hamiltonian strength and 
cooling rates for the direct feedback model. The solid curves are the fidelities of an encoded qubit (three qubit code) with parameters: 
7 = 0.05 Hz, A = 2.5 X k, initial state | Xj/o) = |000). The dashed curves on both figures show the fidelity of one qubit imdergoing 
random bit-flips without error correction (initial state | i//o) = |0)). 



be modeled by the following simple master equation (ME): 

^ = ymxiiii] + mxiii] + mixii])p + x imiis-i] + miiis-])p - mh, pi (3) 

dt 

where 7 is the bit-flip error rate, K is the strength of H, the Hamiltonian evolution which performs the detection and 

correction, and X is the rate of the cooling applied to the ancilla qubits. The super-operator & is defined as in Eq. (2), 
and 5^ = ^{X + iY) = |0) (1 1 is the qubit lowering operator. In the above equation, the ordering of the tensor product 
runs down the circuit (i.e. the first three operators apply to the encoded qubit, and the last two to the ancilla). 

As before, we can evaluate the efficacy of this scheme by numerically solving the ME (3). The results of this are 
summarized by the right subfigure of figure 3, where again, we have used fidelity with the initial state as the figure of 
merit. More details on this scheme and the numerical solutions are in [20]. 

We see from figure 3 that the abihty of the direct feedback scheme to preserve fidelity improves as the Hamiltonian 
strength (k:) and coohng rate (A) are increased. This is to be expected because as these parameters are increased, the 
whole scheme approaches a discrete error correction process (which is the K = X ^ °° limit) operating continually. 
However, it is important to note that we cannot improve performance by simply increasing one parameter; the other 
must also be proportionally increased. 



ASSESSMENT AND CONCLUSION 

The aim of this paper is to show that error correction may be effectively performed using continuous (in time) 
implementations. We have described two continuous implementations that are direct analogies of standard discrete 
implementations, and validated their efficacy through numerical simulations. Two concluding points about the schemes 
presented are: (1) Even though the schemes were evaluated for the bit-flip code, it easy to see how one would extend 
them for any stabihzer code; and (2) Both schemes have a number of free parameters that can be optimized over. We 
expect the performance of such optimized versions to be noticeably better than the results presented here. 

At this point, we can make the following assessment on the advantages and disadvantages of continuous error 
correction implementations: 

• Advantages 

1 . Such models provide a useful framework within which to approach quantum error correction (QEC) from a 
quantum control framework. 



2. Continuous implementations of QEC will be more suited to most quantum computing architectures - at least 
in the short term - because the projective measurements and/or fast unitary gates needed for discrete time 
implementations are not available in most architectures. 

3. Both models can be seen as providing an upper bound for the performance of an error correcting code 
when measurements, gates, and resetting operations are not instantaneous. This is because they simulate the 
situation where such non-instantaneous operations are continuously applied to combat errors - this is an ideal 
scenario because in reality these operations would be interspersed with actual computation operations. 

• Disadvantages 

1. Continuous implementations seem ideal for preserving quantum memory, but are less suited to computation 
because it is not clear that the computation operations can be performed while the error correction procedure 
is on. 

2. In the indirect feedback model, as the QECC being implemented becomes larger, the measurement and 
feedback rates must increase for the scheme to stay effective. This can be a problem practically because 
these rates could be limited by physical restrictions. 

3. In the direct feedback model, the entire set of unitary operations that correspond to the detection and 
correction operations need to be implemented as a Hamiltonian evolution. Depending on the physical system, 
setting up such a Hamiltonian can be a difficult thing to do. 

There are a number of interesting directions in which to continue the investigation of continuous error correction 
implementations. The first of these is to identify good parameter regimes for both of the schemes described above. This 
is important because of the strong dependency of the performance of both schemes on a good choice of parameters. 
Another direction is to investigate how fault tolerance is affected by continuous implementations. Similarly, such 
implementation models could be useful in answering the question: how does the fault tolerance threshold change 
when the operations implementing error correction are not assumed to be instantaneous and arbitrarily accurate? 
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